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FREE FALL AND EVAPORATION OF n-OCTANE DROPLETS IN THE ATMOSPHERE 
AS APPLIED TO THE JETTISONING OF AVIATION GASOLINE AT ALTITUDE 

By Herman H. Lowell 


SUMMARY 

In connection with proposed rapid jettisoning of aviation gasoline 
at altitude, a theoretical investigation was made of the free fall and . 
evaporation of n-octane droplets in the atmosphere. n-Octane was 
selected for study because, on the average, the evaporation characteris- 
tics of that substance approximate those of aviation gasoline. 

Terminal (equilibrium) falling speeds and Reynolds numbers of such 
droplets were obtained for altitudes to 11,000 feet over the droplet 
size range from 6 to 2000 microns. It was concluded that gasoline drop- 
lets larger than 2000 microns (0.2 cm) would be unstable and would there- 
fore occur only rarely. 

The droplet fall data were used in the calculation of evaporation 
rates of falling n-octane spheres; decrease of droplet size and change 
of ambient conditions were taken into account. This procedure was 
applied in the case of a 2000-micron droplet falling at decreasing speed 
from 6000 to about 2000 feet through an NACA standard atmosphere (l5° C 
sea- level temperature). The droplet size decreased to about 200 microns 
during the fall. At that size the falling speed (1.9 ft/sec) was small 
enough and the evaporation rate great enough to ensure subsequent droplet 
disappearance at about the 2000-foot level. 

Additional calculations were made to determine the magnitudes of 
pressure effects, temperature effects, and droplet diameter effects. 

The variation of evaporation rate with pressure for altitudes to 
11,000 feet is very small as compared with the variation with tempera- 
ture and droplet diameter. The evaporation rate increased by a factor 
greater than 100 as air temperature increased from -37° to +30° C. 

It was concluded that use of atomizing devices causing production 
of droplets less than 200 microns in diameter would make possible gaso- 
line jettisoning without ground contamination at ground clearances in 
excess of 250 feet at all temperatures above about -37° C. At the 
highest air temperatures, ground clearances of 1000 feet should suffice 
even when substantial numbers of the largest (0.2 cm) droplets are 
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present, that is, when no atomizing devices are used. At the lowest 
air temperatures, ground clearances as great as thousands of feet will 
not prevent ground contamination if no atomizing devices are used. 


INTRODUCTION 

Suggestions have recently been made that a substantial fraction of 
the fuel load of an airplane should be quickly jettisonable. It has 
been pointed out that fast jettisoning would greatly reduce crash-fire 
hazard and would substantially increase the lift-drag ratio in the case 
of operation under the condition of failure of one or more engines. 

One obvious means of achieving this goal would be the use of dis- 
posable tanks. It is clear, however, that the dropping of fuel-laden 
tanks might not be desirable in many instances. The only practicable 
alternative would appear to be a rapid discharge of the fuel at altitude. 

In general, it is desirable that no significant amount of discharged 
fuel reach any ground installation. The only mechanisms which could con- 
ceivably lead to dispersal at altitude are evaporation and/or break-up 
into minute droplets ("atomization'’) which would remain suspended 
indefinitely. 

The subject of atomization is a very complex one around which an 
extensive literature has accreted. It is assumed here only that in some 
manner, as by the breaking-up of a jet, the fuel has been ejected into 
the atmosphere and is present in the form of droplets of varying sizes. 

The present investigation was conducted at the NACA Lewis laboratory 
with the chief purpose of furnishing preliminary estimates of the frac- 
tions of original masses of droplets of gasoline which will still be 
present after such droplets have fallen through various distances under 
different conditions. It is assumed that n-octane evaporation calcula- 
tions will apply to gasoline. While no estimates are made herein of the 
original distribution of droplet sizes, a brief discussion of the maxi- 
mum stable size and of the motional dynamics of falling gasoline droplets 
is presented in the section MAXIMUM DROPLET SIZE AND DROPLET MOTION. A 
knowledge of droplet speeds is clearly required in the calculation of 
evaporation rates. An explanation is then presented of the rather ele- 
mentary techniques which were employed to obtain the local rates of 
evaporation and the total mass losses. Finally, the results of the cal- 
culations are presented and discussed. 

All symbols used in the report are listed and defined in the appen- 
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MAXIMUM DROPLET SIZE AM) DROPLET MOTION 

Irrespective of the original mechanism of droplet formation^ it is 
intuitively clear that the largest stable droplet during free fall 
through the atmosphere will be that one for which the surface tensile 
forces are just sufficient to overcome the disruptive forces. The latter 
ordinarily consist of the drag shear and normal-pressure forces^ although 
centrifugal forces may also be present (ref. l). Shearing forces may be 
ignored only in situations characterized, by high Reynolds numbers of 
gas-stream flow with respect to the droplet and by rather large droplet 
fluid viscosity (ref. 2). 

Nevertheless, while there are certain contradictions in the litera- 
ture, it is evident that some dimensionless parameter which is essen- 
tially a measure of the ratio of surface forces to deformation-resisting 
forces will determine the maximum stable droplet size under given condi- 
tions. Such a parameter is, for example, the dimensionless Weber number 
defined by the equation 


We 


P^u2r'2 U%' 

oR’ ” a 


( 1 ) 


The effects of tangential forces (shear stresses) cannot he ignored 
at Reynolds numbers below 1000 (ref. 2). It follows that an unknown 
function providing shearing -force information should be added to the 
numerator of equation (l); this fact in no way affects the dependence of 
the resulting modified Weber number on the surface tension. 


The following table, excerpted from reference 2, lists values of 
maximum droplet radius, corresponding free-fall (terminal) droplet speeds 
in normal air, free-fall Reynolds numbers, surface-tension values, and 
values of the ratio (surface tension)/ (critical radius) for several 
substances; 


Liquid 

RS 

c^ 

ft/sec 


a, 

dyne/cm 

a/R', 

dyne/sq cm 

Distilled water 

0.50 

40.02 

16,300 

70.6^ 

141.2 

Carbon tetrachloride 

.18 

30.35 

4,400 

24.5 

136.1 

Methyl salicylate 

.325 

37.07 

9,800 

34.3 

105.5 

Glycerine plus 
20 percent water 

.45 

41.83 

15,300 

63.7 

141.5 


^Probably too low. 
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The approximate constancy of the ratio o/r^ 
average of the four ratios is 131. 


is apparent; the 


In the case of gasoline, the surface tension is about 21 dynes per 
centimeter (ref. 3) . If 131 dynes per square centimeter is accepted as 
a working value, the resulting maximum possible gasoline droplet radius 
is 0.16 centimeter. It is fiirther noted, however, that such a radius is 
rather improbable, as is a water droplet of 0.5-centimeter radius. 
Accordingly, in this report a gasoline droplet diameter of 0.2 centi- 
meter has been arbitrarily selected as representative of the limiting 
size of probably occurring droplets. On the other hand, it should be 
noted that droplets having diameters of the order of 0.3 centimeter will 
occasionally appear; these will be very unstable. 

Equality of drag force and object weight will determine the terminal 
(equilibrium) falling speed of any object. In the case of a (rigid) 
sphere, the following relation applies: 

I gR’3 = Ci)rtR’2 (2) 


This, may be converted into the following: 


32 Pd Pa SR*' 


= Ct 


^Pi 


a 


= Cp Ret^; 


A definition of <?> appears in reference 4 as 

„ ^ 32 Pd Pa 

and is therein presented as a function of Re^ (see fig. 4 of ref. 4)^ 
which may he done since is a unique function of Re^ provided 

the transonic region is not entered. 

Since cp may he calculated for a given radius once the properties 
of the sphere and the amhient medium are known, the terminal Reynolds 
numhers and speeds are easily obtained. 

The density of 115/145 aviation gasoline may he taken, with negli- 
gible error, as 0.700 gram per cubic centimeter (ref. 5). 


(3) 
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The following table exhibits the results of such terminal motion 
calculations for altitudes of 1000, 6000, and 11,000 feet; NACA standard 
atmosphere temperatures and pressures were used: 



Altitude, ft 

1000 

(ta = 55.4° F) 

6000 

(ta = 37.6° F) 

11,000 

(ta = 19.8° F) 


Dd,|i' 

Ret 


Ret 


Ret 

Ut» 

cm 

micron 


ft/sec 


ft/ sec 


ft/sec 

0.1 

2000 

730 

17.41 

700 

19.05 

660 

20.47 

0.03162 

632.4 

90.0 

6.788 

82.5 

7.100 

78.0 

7.648 

0.01 

200 

7.70 

1.837 

7.0 

1.905 

6.40 

1.985 

0.003162 

63.24 

0.3038^ 

0.2894 

0.342^ 

0.2943 

0.3073^ 

0.3013 

0.001 

20 

0.01214 

0.02894 

0.01082 

0.02944 

0.00971 1 

0.03011 

0.0003162 

6.324 

0.0003838 

0.002894 

0.0003421 

0.002943 

0.0003073 

0.003013 


^■For the smaller speeds, Stokes' Law (Cp = 24/Re-t) was used, which is equiva- 
ent to Re-t = <P/24. 


These results are also exhibited in figure 1, in which it is. clearly 
shown that neither terminal Reynolds numbers nor terminal speeds vary 
greatly with altitude within the range covered by the calculations. 

If a droplet is evaporating, the terminal speed and Reynolds number, 
in general, will be constantly decreasing. (An exception to this will 
occur, for example, when- a large, slowly evaporating droplet falls so 
rapidly that the increase in air density with approach to the ground 
more than offsets the diametral change.) The question therefore arises 
as to the validity of the assumption that the instantaneous motion is 
virtually what it would have been under precisely steady-state condi- 
tions, that is, in the absence of evaporative change (or air-density 
change) . The behavior of water droplets under various conditions is 
presented in reference 5, where it is shown that the distances required 
to ensure attainment of equilibrium motion conditions are of the order 
of a few feet for the largest droplets under the most severe accelera- 
tion or deceleration conditions. In none of these cases is the "relaxa- 
tion" distance greater than 40 feet; it is generally less than 25 feet 
for droplets smaller than 0.2 centimeter in diameter. Accordingly, in 
this study it was assiimed that gasoline droplet terminal conditions 
prevailed throughout droplet histories provided the droplets were dis- 
persed sufficiently to ensure absence of aerodynamic interaction among 
them. Since vertical dispersion will occur in any event as a result of 
the large variation of falling speed with diameter, it would appear 
reasonable to neglect such interactions; this procedure will be invalid 
during the initial dispersion period immediately subsequent to the 
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Jettisoning. This initial period will be of doubtful duration, but will 
be brief compared with the total period of droplet evaporation. In any 
case, the assiimption of absence of entrainment phenomena is made through- 
out this study. 


EVAPORATION CALCULATIONS 

The rate of evaporation (mass loss per unit surface area per unit 
time) is equal to the product of two variables, namely, the mass-transfer 
coefficient and the difference in vapor pressure between droplet and 
ambient atmosphere. The droplet vapor pressure is a unique function of 
surface temperature, while the ambient vapor pressure may be taken as 
essentially zero in the atmosphere (except possibly during the initial 
dispersion period). 

The droplet temperature will depend, as in the case of the common 
wet-bulb psychrometer, on the location on the droplet temperature scale 
of the balance point between sensible heat transferred to the droplet by 
conduction and convection from the atmosphere and heat removed from the 
droplet, principally by the evaporating fluid. 

The heat- and mass -transfer coefficients depend on the droplet size 
and on the instantaneous relative aii*-mass flow rate, that is, the drop- 
let falling speed and local air density. Since periods of free fall in 
the atmosphere are long compared with the rates at which successive 
equilibriimi states are attained under various conditions, it is per- 
missible to assume that the droplet is, at every point, essentially in 
psychrometric equilibrium with the local environment. 

Most of the calculations of local mass-loss rates were carried out 
on a step basis. In outline form, the procedure was as follows: 

The initial falling speed and Reynolds number having been deter- 
mined for the starting altitude (6000 feet; NACA standard atmosphere) 
and assumed droplet size (0.2 cm), the initial heat- and mass-transfer 
coefficients were computed. In turn, the droplet temperature, vapor 
pressure, and mass-loss rate were computed. The necessary relations 
are subsequently given. 

The duration of the first step was then taken as that period during 
which a fixed, predetermined fractional mass loss (0.10875) occurred. 

The fractional loss was based on the mass existing at the end of the 
preceding step, rather than on the original mass. Therefore, each 
successive mass was 0.89125 of the preceding mass. 
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The diameters and masses of the successive steps were thus estab- 
lished in advance. The curves of figure 1 reveal that falling speeds 
and Reynolds numbers are nearly independent of altitude within the alti- 
tude range studied. It was therefore possible to establish beforehand 
a list of terminal speeds and Reynolds numbers corresponding to the 
successive drop sizes for the starting altitude. The diameters used in 
the calculations were values more nearly representative of average con- 
ditions existing over the successive intervals than either the respective 
initial or final values would have been. The mean mass was taken as 
0.945625 of the starting mass at each step (0.945625 = 1 - 0.10875/2). 

The mean diameter was then taken as 0.945625^ or 0.981534 of the start- 
ing diameter. (If the mean diameter is defined by the relation 

- 2 _2 -2 

= 0.5(D(i^i + ), where Dd.. i is the mean diameter over 

the i"^^ step and and D^^i+l are the respective initial and final 

diameters, the mean diameter so obtained is also about 0.981.) 

Once the period of fall and rnesn speed had been determined, the loss 
of altitude for the first step was then available. A more nearly exact 
mean temperature and pressure for that step were than obtained from 
reference 6j the step calculations were then repeated. After the condi- 
tions at the end of the first step had been calculated, the distance of 
fall for the second step was estimated, the average temperature and 
pressure were obtained, and the procedure of the first step repeated. 

The entire cycle was reiterated until the desired droplet size was 
reached. 

The relations used were the following; 

Heat-balance or psychrometric equation: 


^a - "td = (Pd - Pa) (4) 

It has heen shown hy many investigators (ref. 7 , for example) that 
this relation is valid provided that no substantial amount of sensible 
heat is transferred by the liquid vapor. This will be true, in general, 
when both the mass- and heat-transfer rates are low. Further, if this 
is the case, the mass- and heat-transfer rates will be essentially 
independent . 

Under such conditions, mass-transfer and heat-transfer Nusselt 
numbers may be defined (refs. 4 and 7) as follows: 
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Num s 


NUhH 


P 




( 5 ) 




J 


For most pairs of gases, the diffusion coefficient 8 varies with 
temperature and pressure in the following manner (ref. 8): 


P = fo (^) (6) 

The question arises as to the constancy of Pq of* other liquid or 

gas properties during an evaporation process. 

In the case of gasoline, a complex mixture of many compounds, a 
fractional distillation process, in effect, takes place as evaporation 
proceeds. Accordingly, such properties as the air-gasoline vapor 
diffusion coefficient, heat of vaporization, and so forth will vary 
markedly with time. 

A rather radical solution of the difficulty consists in replacing, 
on paper, the gasoline droplet by a droplet of liquid which, on the 
average, exhibits roughly the same evaporative behavior as gasoline. 

The term "evaporative behavior" refers here only to the rate of loss of 
mass and not, for example, to the properties of any of the successive 
fractions evaporated. The selection of n-octane as a representative 
single material should yield evaporation data reasonably similar to 
those for gasoline. The correspondence will not be close during the 
initial evaporation period. At that time, the higher fractions, which 
have vapor pressures a number of times as high as that of n-octane, 
will evaporate at a high rate. Further, after most of the original 
gasoline has been vaporized, the converse will be true; the evaporation 
rate of the residual liquid will be lower than that of n-octane. 

Nevertheless, the results of calculations based on the properties 
of n-octane are believed to be a;pproxiinately correct for a mixture 
including substantial quantities of n-octane, isooctane, heptanes, and 
other paraffins. Accordingly, the value of Pq cited in reference 8 
for n-octane. was used, namely, 0.0505 square centimeter per second. 
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The following additional constants were used in the calculations; 
the value of the Schmidt number Sc,! discussed subsequently: 


Constant 

Magnitude 

Units 

Source 

Thermal conductivity 
of air, kh,o 

2422 

erg/ cm 

Reference 9 

Molecular weight of 
air, Ma 

28.966 

gm/mole 

Reference 10 

Molecular weight of 
n-octane, M^ 

114.224 

gm/mole 

Reference 11 

Universal gas 
constant, R 

8.3144X10'^ 

erg/°K 

Reference 11 

Schmidt number for 
n-octane (in air), 

Scd,o 

2.618 



Heat of vaporization 
of n-octane, X^^q 

3.734X109 

erg/gm 

Reference 12 data on vapor 
pressure and extended plot 
of heat of vaporization 
data in reference 13. 


The following substitutions are now made in equation (4): (l) hjjj 

and hjj are replaced by their equivalents as obtained by solving equa- 
tion (5); (2) p is replaced in the expression for hm by relation (6); 
(3) kh is replaced in the expression for h^ by = k^ e(Ta,/N.T. 
and (4) is replaced in the final expression by X^ = X,^^g(X,!/X,!^g) . 

The resulting expression is the following: 



In equation (7), the temperature difference is given in degrees 
centigrade, while the vapor pressure difference is given in dynes per 
square centimeter. The vapor pressure is generally given in millimeters 
of mercury in the literatiire; if p* is used to represent a vapor pres- 
sure so expressed, equation (?) becomes the working form of the psychro- 
metric equation for n-octane under moderate ten^ierature and pressure 
conditions: 
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'Nu_ 

ta - td = 0.522 ( 


h 


N.P. 

Po 


N.T. 


0.15 




(Pd - pD (8) 


This states that, at low relative flow rates, the temperature drop 
at normal pressure of an p-octane sphere is roughly half the difference 
in ^por pressxire existing between the sphere surface and the atmos- 
phe^, when the latter is expressed in millimeters of mercury. In most 
cases, p* will be virtually zero. The reason for the restriction of 
the observation to low relative flow rates lies in the variation of the 
ratio Nti^/Nujj, which can only be experimentally determined. 


The following expressions have been obtained from data cited in 
reference 4 by a trial-and-error fitting process; it is not contended 
that these are the "best" expressions possible, but only that they fit 
the data fairly well for moderate (^1000 ) Reynolds numbers: 


1 ^ 
= 2 + 0.39 Sc^ Re°*^® 

1 

Wuj^ = 2 + 0.246 Pr^ 


( 9 ) 


These expressions are represented graphically in figure 2, while 
their ratio NvIj^/Nur is presented in figure 3. Fourth nonzero figures 
are not of significance here; third significant , figures may be erroneous 
by several units. The values plotted, however. Were those used in the 
calculations. 


The value 0.73 was adopted as the Prandtl number for pure air. The 
Schmidt number is given in reference 7 as 2.58 at 25° C; this value 
agrees with a value calculated on the basis of a diffusion coefficient 
of 0.060 at 25° C, which is the value cited in the same reference. In 
turn, the value 8 = 0.060 agrees well with the value obtained by using 
equation (6) and Pq, namely, 0.06010. If the latter is used, the Schmidt 
number at 25° C is found to be 2.57. Finally, the Schmidt number, which 
varies very slowly with temperature, was 2.618 at 0° C for n-octane, 
when Pq = 0.0505 and the relation Sc^ s Pa/PaP are used.” 

Vapor pressure was plotted against heat of evaporation (fig. 4) to 
facilitate trial-and-error solutions of equation (8) . In the original 
figure, the temperature could be read to a tenth of a degree and the 
vapor pressure to at least 1 percent. 

As a check on the validity of relations (8) and (9), they were used 
to calculate wet-bulb temperatures of n-octane spheres at three air 
temperatures. The results for a Reynolds number of 1000 are compared in 
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figure 5 with experimental data from reference 14 obtained at about the 
same Reynolds niimber. The agreement is to some extent fortuitous, but 
it is clear that, at least under roughly normal ambient conditions, the 
theoretical heat- and mass-transfer rates calculated from the relations 
given here will be substantially correct (for n-octane). 

The mean fractional evaporation rate at the i"^^ step (designated 
g/(s)(sec)) was calculated as follows; note the assumption in equa- 
tion (lO) that the partial n-octane vapor pressure in the atmosphere 
is negligible: 


i 1 ■%, i Pd, i 


= jt D, .2 h . p, . 
d, 1 m, 1 -^d, 1 


(10) 


It therefore follows that 


Ai = 


jt D , . h . Pj . 

d, 1 m, 1 d, 1 


(n/6) Pd Dd,i" 


= ^ Vi 

pd^^T 

The following substitutions and changes are now made: 


“ 1333.2 
“ 0 • 7 00 




Md 

% ^ HTa 


I^] 

VN.T.y 



Md 

RTa 


2 
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or 


T— — Po /Ta, iWn.P. 

5^ (H.T.)Ea 


xio^ 


( 11 ) 


The fractional mass loss rate then becomes 


= I 6X1355.2X10^0.0505 \ / N.P." ^ 

\ 0.700X273.2X0.7279X10®/ 5^^^ / V^a,i> 




or 



Nu . 

2.9oxio2 : 

\x, i 




( 12 ) 


The number of seconds required for a fractional loss of 0.10875 
of the mass present at the beginning of a step is then 


Ai0 


0.10875 


(13) 


The distance traveled during the i^^ step is the product of A±9 
and Ut^iJ total distance fallen and the existing altitude are then 
obtainable in an obvious manner. 


' RESULTS AND DISCUSSION 

The original intention had been to have complete sets of calcula- 
tions performed for the fall of droplets of several sizes from several 
altitudes at several temperatures. The achievement of a reasonable 
degree of accuracy in the calculations required, however, computations 
somewhat lengthier than at first thought necessary, and the program was 
therefore abridged. 

Calculations for the case of a droplet having an initial diameter 
of 2000 microns and falling through a standard atmosphere from 6000 feet 
have been carried sufficiently far at this time to make possible pre- 
dictions concerning the release of gasoline from the altitude in question 
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The mass of the droplet as a function of altitude is shown in figure 6; 
the respective diameters are given at relative droplet mass values of 
1^ 0.1, 0.01, and 0.001. Only the results of calculations at each fifth 
step have heen plotted. 

It is clear that large droplets of n-octane (and, therefore, 
probably of gasoline) will fall large distances (at least 4000 ft) under 
normal atmospheric conditions before their sizes are reduced to values 
such that the droplets will have essentially come to rest. It is to be 
noted (fig. 1) that gasoline droplets (or n-octane droplets having about 
the same density) of 200 microns diameter will fall freely (at 6000 ft) 
at a speed of 1.9 feet per second. It is safe to assume that such drop- 
lets will not reach the ground unless air temperatures are extremely 
low; it is estimated that a 200-micron droplet will lose 11 percent of 
its mass while falling 54 feet at -37^ C (= -35^ F). 

Figure 6 may be used to determine residual masses of smaller drop- 
lets falling from the altitudes indicated along the abscissa scale. For 
example, at 4696 feet the residual mass of the original droplet is 
0.5623 of the original mass (which was 2.932X10"^ g); at 2972 feet the 
residual mass of the original droplet is 0.1000 of the original mass. 
Further, the respective diameters at the two altitudes mentioned are 
1651 microns and 928 microns. It is clear, considering the assumptions 
which underlie the calculations, that the residual mass of a 1651-micron 
droplet falling from 4696 feet to- 2972 feet (in normal atmosphere) , will ' 
be 0.1000/0.5623 of the original mass, or 0.1778 of the original mass. 

The first of the following tables exhibits in abridged form the 
respective masses, diameters, and so forth, for the set of calculations 
in question; the second table exhibits the masses, diameters, and so 
forth, for a small number of calculations which were made in lieu of . 
additional complete sets of calculations to detennine the effects of 
initial diameter variation at a fixed altitude, temperature variation 
at. a fixed pressure altitude, and altitude variation at a fixed tem- 
perature: 
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Step 

num- 

ber 

Mi/Ml 

and 

microns 



°C 

^d^ 

Ai> 

g/(g)(sec) 

sec 

Alti- 

tude^ 

ft 

Ret,i 

ft/sec 


1.0000 

2000 







6000 

1 


1963 

670 

18.9 

3.43 

0.67 

0.00653 

16.65 



.5623 

1651 







4696 

6 


1620 

478 

16.2 

5.89 

2.79 

.00901 

12.07 



.3162 

1363 







3869 

11 


1337 

340 

13.9 

7.46 

4.13 

.01194 

9.11 



.1778 

1125 







3330 

16 


1104 

237 

11.8 

8.49 

4.99 

.01530 

7.11 



.10000 

928.3 







2972 

21 


911.1 

167 

10.0 

9.17 

5.54 

.01962 

5.54 



.05623 

766.2 







2734 

26 


752.1 

115 

8.41 ! 

9.63 

5.94 

.0247 

4.41 



.03162 

632.5 







2576 

31 


620.8 

79.8 

7.00 

9.92 

6.20 

.0314 

3.47 



.01778 

522.0 







2472 

36 


512.4 

54.0 

5.81 

10 . 12 

6.40 

.0394 

2.76 



.01000 

430.9 







2405 

41 


422,9 

•36.4 

4,73 

10.25 

' 6.54 

.0503 

2.16 



.00562 

355.7 







2362 

46 


349.1 

24.3 

3.82 

10.33 

6.68 

.0640 

1.70 



.00316 

' 293.6 







2335 

51 


: 288.1 

16.2 

3.05 

10,38 

6.79 

.0826 

1.32 



.00178 

242.3 







2319 

56 


247.1 

11.4 

2.52 


— 





Note: Values given on line above that on which step number is entered are 

those applying to conditions at the beginning of step. Values on step 
line are average values. 


Diameter effects 

Altitude, 


Ret 

Ut, 

^a> 


A, 

A8, 

AH, 

ft 

microns 


ft/sec 


°c 

g/(g) (sec) 

sec 

ft 

6000 

1963 

670 

18.9 

3.11 

0.39 

0.00642 

16.94 

320 


911.1 

167 

10.0 


.28 

.01500 

7.25 

72.5 


422.9 

36.4 

4.73 


.26 

.0376 

2.89 

13.68 


196.3 

6.62 

1.84 


.54 

.1044 

1.042 

1.917 


91.11 

.940 

.554 


.81 

.345 

.315 

.1746 


42.29 

.107 

.137 


.99 

1.402 

.0776 

.0108 

Temperature effects 


6000 

1963 

670 

18.9 

-37.00 

-37.12 

0.000226 

481 

9100 





-12.00 

-13.10 

.00246 

44.3 

837 





3.11 

.39 

.00642 

16.94 

320 





13.00 

8.30 

.01151 

9.45 

178.6 





30.00 

21.0 

.0234 

4.66 

88.0 

Altitude effects 

1000 

1963 

710- 

17.1 

13.00 

8.99 

0.01011 

10.76 

184.0 

6000 


670 

18.9 


8.30 

.01151 

9.45 

178.6 

11,000 


638 

20.0 


7.74 

.01251 

8.69 

173.9 





NACA RM E52L23a 


15 


The droplet diameters used were average values corresponding to 
selected steps of the step-calculations. On the other hand, the pres- 
sures used were those of the NACA standard atmosphere at the indicated 
altitudes rather than the mean values over an initial step. 

In the previous table, selected results are given for the case of 
a 2000-micron droplet falling through an atmosphere, the characteristics 
of which change with altitude. The "Diameter effects" group of the sec- 
ond table exhibits differences, to the contrary, caused solely by diam- 
eter changes. A detailed comparison with the original step calculations 
for the extended-fall case indicates that the higher evaporation rates 
of that situation may easily be accounted for on the basis of the higher 
air temperatures existing at the lower altitudes. For example, in the 
step calculations (previous table), the number of seconds required for 
the standard 11-percent mass loss to occur at step number 21 is 5.54 
(at a mean altitude of 2940 ft) . In the present case, the nvimber of sec- 
onds required for a droplet of the same diameter (911.1 microns) to lose 
that much mass Is 7.25. The respective droplet temperatures are 5.54° C 
(a coincidental repetition of the same number) and 0.28° C. The wet-bulb 
depression is roughly the same for both cases, namely, some 5° or 6° C, 
while the ratio of the two evaporation rates is roughly proportional to 
the ratio of the two vapor pressures. Such reasoning will facilitate the 
estimation of evaporation rates under varying conditions; temperature and 
diameter effects may be combined in estimating loss rate for a particular 
droplet size at a particular temperature. 

The very great variation of evaporation rate with temperature is 
indicated in the second group of this table. The distance of fall at 
-37° C is so large that the assumption of constancy of temperature is, of 
course, not tenable; it is clear, however, that under severe conditions 
reliance may not be placed on evaporation as a dispersing mechanism. On 
the other hand, the loss rates at high temperatures are such that moder- 
ately low-altitude jettisoning even without the use of finely atomized 
sprays would be feasible, since the largest droplets (2000 microns) can 
drop, at most, perhaps 1000 feet before essentially coming to rest. 

(This estimate is based in part on the 4000-ft drop previously found for 
the case of a 2000-micron droplet falling through an NACA atmosphere.) 

The evaporation figures of the "Altitude effects" group show quite 
clearly the rather small effects of pressure differences over the indi- 
cated altitude range. 


CONCLUSIONS 

A theoretical investigation was made of the free fall and evapora- 
tion of n-octane droplets in the atmosphere. On the basis of the analy- 
sis and numerical results the following conclusions were drawn; these 
conclusions will be approximately applicable to the fall of aviation 
gasoline : 
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1. As the droplets evaporate, their masses, diameters, and there- 
fore their falling speeds decrease. Motion and evaporation analyses are 
therefore most easily made on a step basis, with changes of pressure and 
temperature with altitude taken into account; droplet characteristics 
and environmental conditions are supposed constant during each step. 

2. Droplets having diameters greater than 0.3 centimeter will either 
not be formed or, having been formed, will break up immediately into 
smaller droplets. The largest stable droplets will have diameters of 
the order of 0.2 centimeter. 

3. Temperature effects (at a given droplet diameter) on evaporation 
rates, and therefore on total distances of fall, are far more important 
than pressure or motion effects. The rate of evaporation increases by a 
factor of more than 100 over the air-temperature range -37° C (s -35° F) 
to +30° C (86° f). 

4. A droplet having a diameter of 200 microns (0.02 cm) or less 
will fall at a speed of less than 2 feet per second at altitudes below 
11,000 feet. The rate of evaporation, and therefore of decrease of size 
and falling speed, of any such droplet will be great enough even at air 
temperatures as low as -37° C to ensure that the droplet will not reach 
ground level when ground clearance exceeds about 250 feet. It follows 
that ground contamination will not occur, irrespective of air tempera- 
ture, if nearly all of the jettisoned fuel has been atomized, that is, 
if the fuel is present in the form of droplets less than 200 microns in 
diameter at ground clearances in excess of 250 feet. 

5. Large droplets will fall long distances through air under normal 
or near-normal conditions. For example, in an NACA standard atmosphere 
(15° C sea-level temperature), a droplet having a diameter of 2000 microns 
(0.2 cm) at 6000 feet will fall while evaporating at decreasing speed to 
about 2000 feet. At this level it would have a diameter of less than 

200 microns and would remain approximately at this level while evaporating 
completely. 

6. At the highest air temperatirres (about 30° C), ground clearances 
of 1000 feet will probably be sufficient to avoid ground contamination by 
falling droplets, irrespective of the original distribution of droplet 
sizes. At the highest air temperatures, therefore, the use of finely 
atomizing spray devices will be unnecessary. 

7. At the lowest air temperatures (below about -20° C), only the use 
of finely atomizing spray devices will prevent ground contamination even 
if ground clearances are of the order of thousands of feet. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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APPENDIX - SYMBOLS 

A surface area^ sq cm 

Cj) drag coefficient for a sphere 

D diameter^ cm 

g gravitational constant, dyne/g 

H altitude, cm (except where explicitly expressed in ft) 
hji heat-transfer coefficient, erg/(sq cm) (sec) (^C) 
hm mass-transfer coefficient, (g/(sec) (sq cm) )/(dyne/sq cm) 
kjj thermal conductivity of ambient gas, erg/ (cm) (sec) (°C) 

M molecular weight, g/mole 

Nu Nusselt number 

N.P. 1 atmosphere, dyne/sq cm (1,013,250) 

N.T. reference temperature, °K (273.2) 

P ambient pressure, dyne/sq cm 

Pr Prandtl number 


Pa 

Pd 

Qm 

R 

Rd 

Re 

r 

Sc 

T 


partial pressure of evaporating- liquid vapor in ambient gas, 
dyne/sq cm 

saturation vapor pressure of evaporating liquid at sphere tempera- 
ture, dyne/sq cm 

rate of mass loss from sphere, g/sec 

universal gas constant, erg/ (°K) (mole) 

gas constant of evaporating liquid, erg/(°K)(g) 

relative droplet -air Reynolds number 

droplet radius , cm 

Schmidt nvuaber, Pa/pg,? 

absolute temperature, °K 
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t temperature, °C (except where explicitly expressed in °F) 

U relative air-droplet speed, cm/sec (except where explicitly 
expressed in ft/sec) 

We Weber number defined by equation (l) 

a surface tension, dyne/cm 

3 diffusion coefficient, sq cm/sec 

A fractional mass-loss rate, g/(g)(sec) 

A heat of vaporization, erg/g 

H viscosity, poise 

p density, g/cu cm 

( 

<p free-fall parameter defined by relation following equation (3) 
Subscripts; 
a ambient gas 

I c critical value (referring to value above or beyond which droplet 
tends to break up) 

d evaporating liquid 

h heat transfer 

i step number designation 

m mass transfer 

o reference conditions. (N.P., W.T.) (Note that these are not identic 
cal with NACA standard atmosphere sea-level conditions. In the 
latter case, the air temperature is taken as 15° C.) 

t terminal (equilibrim) state 

p dimensions microns instead of cm 

The bar* (— ) denotes mean value over a step (calculation interval) . 
Superscripts: 

* dimensions mm Hg instead of dyne/sq cm 
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Droplet diameter, microns 

Figure 1. - Terminal motion characteristics of gasoline droplets. 
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Figure 3. - Ratio of mass- and heat-transfer Nusaelt numbers for n- octane sphere evaporating Into normal air as function of Reynolds number. 
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Figure 5 , - Psychrometrlc behavior of n-octane spheres . Reynolds number equal approximately to 1000 . 
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Figure 6. - Mass loss by evaporation of 2000-mlcron n-octane droplet falling through NACA standard 
atmosphere . 
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